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1. Introduction 

In Ref. [1] , we presented a formalism for a mathematical representation of a parton shower 
that incorporates interference in both spin and color. In this paper, we analyze this formal- 
ism in the approximation that we average over parton spins at each step and keep only the 
leading contributions in an expansion in powers of where A'^c = 3 is the number of 

colors.^ Our interest is to elucidate the structure of the full shower formulation of Ref. [1] 
by examining what happens when the spin-averaged and leading color approximations are 
imposed. We also anticipate that the approximate shower may be of use in implementing 
successively better approximations to the full shower including spin and color. 

^More precisely, we average over the spins of incoming partons at each step and sum over the spins of 
the outgoing partons. 
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Our main focus is on the splitting functions that would be used to generate the shower 
in the spin averaged approximation (which is a customary approximation in current parton 
shower event generators). In our formalism, there are two sorts of splitting functions. The 
direct splitting functions correspond to the squared amplitude for a parton / to split into 
daughter partons that, in our notation, carry labels / and m + 1, where m + 1 is the total 
number of final state partons after the splitting. In this paper, we use the spin dependent 
splitting functions from Ref. [1] and simply average over the spins of the mother parton and 
sum over the spins of the daughter partons. We analyze some of the important properties 
of these functions. We also need interference splitting functions. These correspond to the 
interference between the amplitude for a parton / to split into partons with labels I and 
m + 1 and the amplitude for another parton k to split into partons with labels k and 
m + 1. These functions generate leading singularities when parton m + 1 is a soft gluon. 
We improve the specifications of Ref. [1] for this by defining a useful form for certain weight 
functions Ai^ and Ai^i that were assigned the default values 1/2 in Ref. [1]. We will see that 
with the improved form for Aij, the total splitting probabilities acquire useful properties 
in the soft gluon limit. 

We will see that when we make the spin-averaged and leading color approximations, the 
parton shower formalism of Ref. [1] amounts to something quite similar to standard parton 
shower event generators. One significant point in common is that the splitting functions are 
positive. One difference with some standard event generators is that an angular ordering 
approximation is not needed because the coherence effects that lead to angular ordering are 
built into the formalism from the beginning, both for initial state and final state splittings. 
This coherence feature is a natural consequence of a dipole based shower, as in the final 
state showers of Ariadne [2] and the kx option of Pythia [3] or the showers [4, 5] based 
on the Catani-Seymour dipole splitting formalism [6]. Additionally, our formalism differs 
from others in its splitting functions and its momentum mappings. 

2. Direct spin-averaged splitting functions 

We begin with the splitting functions that correspond to the amplitude for a parton to 
split times the complex conjugate amplitude for that same parton to split. We follow 
the notation of Ref. [1]. Before the splitting, there are partons that carry the labels 
{a,b, 1, . . . ,m}, where a and b are the labels of the initial state partons. The momenta 
and flavors of these partons are denoted by {p, f}m = {Psl, fa] ■ ■ ■ ]Pm, fm}- The flavors 
are {g,u, u,d, . . . }, with the initial state flavors /a and /b denoting the flavors coming out 
of the hard interaction and thus the opposite of the flavors entering the hard interaction. 
We let I be the label of the parton that splits. After the splitting, there are m + 1 final 
state partons. The momenta and flavors of the partons are {p, f}m+i- We use the label / 
for one of the daughter partons and the label m + 1 for the other daughter parton.^ The 
partons that do not split keep their labels. However, they donate some of their momenta 
to the daughter partons so that the daughter partons can be on shell. Thus pi ^ pi in 

^For a final state q ^ qg splitting, we use m + 1 for the label of the gluon. For a final state g ^ qq 
splitting, we use m + 1 for the label of the q. 
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general for a spectator parton. The momenta and flavors after the sphtting, {p, f}m+i, are 
determined by the momenta and flavors before the sphtting, {p, f}m and variables {Cp, Cf} 
that describe the splitting.^ A certain mapping 

{p, fU+i = Ri{{p, fU, {Cp, Cf}) (2.1) 

deflned in Ref. [1] gives the relation. 




Figure 1: Illustration Eq. (2.2) for a qqg final state splitting. The small filled circle represents 
the splitting amplitude vi. The mother parton has momentum pi + Pm+i, but in the amplitude 
|Af({p, /}m)), this ofF-shell momentum is approximated as an on-shell momentum pi. 

The splitting functions in Ref. [1] are based on spin dependent splitting amplitudes vi. 
One starts with the amplitude |M({p, /}m)) to have m partons. The amplitude is a vector 
in spin(8)color space. After splitting parton /, we have a new amplitude \Mi{{p, f}m+i)) of 
the form illustrated in Fig. 1 

\Mi{{p, fU+l)) = tjifi ^fi + fm+l) ^/({P, fU+l) \M{{p, f}m)) . (2.2) 

Here tj is an operator on the color space that simply inserts the daughter partons with 
the correct color structure. The factor V^^ is a function the momenta and flavors and 
is an operator on the spin space. It leaves the spins of the partons other than parton 
/ undisturbed and multiplies by a function vi that depends on the mother spin and the 
daughter spins: 

{{s}m+l\yi{{p,f}m+l)\{s}m)={ dg^^sA Vl{{p, f}m+l, Sm+1, k, Si) . (2.3) 

\j(^{l,m+l} J 

Thus the splitting is deflned by the splitting amplitudes vi, which are derived from the 
QCD vertices. 

^When a gluon splits, (^f determines wfiether the daugliters are a (g, g) pair, a (u, u) pair, etc. In Ref. [1], 
we defined the sphtting variables (p in a rather abstract way, but one could imagine using for the 
virtuality of the daughter parton pair, a momentum fraction, and an azimuthal angle. 
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We can illustrate this for the case of a final state q ^ q + g splitting, for which we 
define 



Vl{{p, f}m+l,Sm+l, Sl, Si) = 

U{pi, h)l^\Pi + fm+i + m{fi)\rliU{pu si) (2.4) 



V4vras e^(pm+i, Sm+i; 0)* 



[{Pi +Pm+iY - m'^ifi)] '^Pi-n-i 



There are spinors for the initial and final quarks. There is a polarization vector for the 
daughter gluon, defined in timelike axial gauge so that Pm+i ■ e = Q ■ e = 0. Here Q is the 
total momentum of the final state partons, which is the same before and after the splitting. 
There is a vertex for the qqg interaction. There is a propagator for the off shell quark 
with momentum pi + pm^i. So far, this is exact. Finally, there is an approximation that 
applies when the splitting is nearly collinear or soft. We approximate pi + Pm+i by pi in 
the hard interaction and insert a projection ifLj2pi ■ rii onto the "good" components of the 
Dirac field. This projection uses a lightlike vector ni that lies in the plane of Q and pi, 

ni = Q ^ pi . (2.5) 

Q-Pi + ViQ-Pi)'-Q'^''ifi) 




Figure 2: Illustration of how vi times Uj* appears in the calculation of the approximate matrix 
element |M;({p, /}„j+i)) times its complex conjugate, (M;({p, If we average over spins, 
we need to multiply \M{{p,f}m)) times (M({p,/}„)| by Wu, Eq. (2.8). 

With one exception, the direct splitting functions in Ref. [1] are products of a splitting 
amplitude, vi, times a complex conjugate splitting amplitude, u^*, 

Vl{{pJ}7n+l,Sm+l,Sl,Sl)vi{{p,f}.m+l,s'^+l,s'i,s'i)* . (2.6) 

The calculation of \Mi{{p, f}m+i)) times (^Mi{{p, f}m+i)\ using vi x wj* is illustrated in 
Fig. 2. In this calculation, in general, we have to keep track of two spin indices, s and s' for 
each parton in order to describe quantum interference in the spin space. However, in this 
paper we make an approximation. We set s' = s for each parton, sum over the daughter 
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parton spins and average over the mother parton spins. Thus we use a sphtting function^ 



1 



\vi{{p, /}m+l, k, Sl)Y 



Sl,Sm + l,Sl 



for any flavor combination ahowed with our conventions for assigning the labels / and m + 1 
except for a final state g ^ g + g splitting, for which we do something slightly different 
because the two gluons are identical. We make manifest the definition of which flavor 
combinations are allowed by defining 



Sl{{f}m+l) 



r 1/2, /e{i, 
1, ^e{i, 

0, /G{1, 

1, /e{i, 

0, /E{1, 
I 1, ZG{a,b} 



:'m}, fi = fm+l = g 

fl / gjrn+l = g 
fl = gjm+l / g 

^m}, fl = qJm+1 = q 
fl = qJm+1 = q 



(2.7) 



This is 1 for the allowed combinations, otherwise, with a statistical factor 1/2 for a final 
state g ^ g + g splitting. The complete definition of Wn is then 



Wu 



Sl{{f}m+l) ^ ^ l\vi{{p,f}m+l,Sm+l,Sl,Sl)\'^ 

Sl,Sm + l,Sl 

+ e{le{l,...,m}Ji = f,n = g) 

X \v2,l{{pJ}ni+l, Sm+1, Si, Sl)]"^ - \v3,l{{pJ}m+l, Sra+1, h, Sl)f 



The second term applies for a final state g ^ g + g splitting and is arranged to keep the 
total splitting probability the same but associate the leading soft gluon singularity with 
gluon m + 1 rather than gluon I. The functions V2,i and v^^i are defined in Sec. 2.3. 

The form of the splitting amplitude vi depends on the type of partons that are involved. 
However, there is a common result in the limit Pm+i whenever parton m + 1 is a gluon. 
In this limit, vi is given by the eikonal approximation. 



, .eikonal 



Vl 



{{P,f}m+l,Sm+l,Sl,Sl) = y/47TQ^6si,si 



^{Pm+l,Sm+i;Q)*-Pl 
Pm+l -Pl 



(2.9) 



The soft gluon limit of Wu is then 



W 



eikonal 
II 



Pl ■ D{pm+V,Q) -pi 

{Pm+I-Plf 



(2.10) 



Here D^^ is the sum over Sm+i of e^e^. 



Pm+lQ^ + Q^Pm+l Q'^Pm+lPm+l 



D^u{Pm+l,Q) = -QiJLU + 



Pm+l ■ Q 



{Pm+1 ■ QY 



(2.11) 



*The function Wu here is the same as wu in Ref. [1]. 
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The function Wu and its approximate form W^^"^^^ give the dependence of the sphtting 
operator on momentum and spin for a given set of parton flavors. The partons also carry 
color. In Ref. [1] there is a separate factor that gives the color dependence. This factor is 
an operator on the color space that we can call tj ^ti, where tj is the operator in Eq. (2.2), 
which inserts the proper color matrix into the amplitude, and inserts the proper color 
matrix into the complex conjugate amplitude.^ So far, we do not make any approximations 
with respect to color. In Sec. 5, we will make the approximation of keeping only the leading 
color conributions. 

We now turn to a more detailed discussion of Wu for particular cases. 

2.1 Final state + g splitting 

Let us look at Wu for a final state g — > g + g splitting, 

47roi 1 

= i2PrP^,.r ""^^^^-^^ (2.12) 



X -Tr 
4 



Here m = m{fi) is the quark mass, the lightlike vector ni is given by Eq. (2.5), and D^^, is 
given by Eq. (2.11). It will be convenient to examine the dimensionless function 

F ^ PlA^ Wu . (2.13) 
Anas 

The limiting behavior of F as the gluon m + 1 becomes soft, Pm+i 0, is simple. Then 
the eikonal approximation applies and we obtain from Eq. (2.10) 

F..onal = ^"^;^-""^^-^' . (2.14) 
Pi ■ Pm+1 

The full behavior of F is more complicated, 

F=[l + h{y, ai, bi)] Feikonal + P^^^±^ . (2.15) 

Pi ■ ni 



Here 



where 



h{y,ai,bi) = — — h , , ^-1, (2.16 

1+ri Xri{l+ri) 



"^Pi ■ Pm+l 
2pi-Q 



ai 

bi 
ri 
A 



2pi-Q 

(2.17) 

2pi-Q 

\/l - 4a/ 6; , 

^{l + yY -Aai{y + bi) 

n 



^In Ref. [1], we write t\{fi fi +g) for the operator that we here call just t^^ and we denote the operator 
tj(g)ti hyg{l,l). 
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The eikonal approximation to F will turn out to be significant in our analysis when 
we incorporate the effect of soft-gluon interference graphs. We will find that it is of some 
importance for the numerical good behavior of the splitting functions including interference 
that 

F - Feikonal > . (2.18) 

To see that this property holds we note first that Pm+i ■ ni/pi ■ ni is non-negative. Remark- 
ably, h{y,ai,bi) > also. To prove this, we first note that 

h{y,ai,0) = — , (2.19) 

so that h{y,ai,0) > 0. Then we show that h{y,ai,bi) — h{y,ai,0) > by simply making 
plots of this function. This establishes the positivity property Eq. (2.18). 

We now examine F further under the assumption that m = 0. We write F as a function 
of the dimensionless virtuality variable y, and a momentum fraction^ 



^ _ Pm+i ■ ni 
{pm+i + Pi) ■ ni 

It is also convenient to use an auxiliary momentum fraction variable 



(2.20) 



P-^' • ^ ^ z+ , - , (2.21) 



(p^+i+k)-0 1 + y (l + y)(l + y + A) ' 
where, for m = 0, A = y^(1 -|- y)^ — 4a;y. Using these variables, 

\ — X '^CLiy 

-feikonal — 2 ■ , (2.22) 

X X^(l -I- y)^ 



and 

F 



{X-l + yf+4y 
4A 



i^cikonai + ^z[l + y + X] . (2.23) 



As y ^ 0, F must turn into the Altarelli-Parisi function for this splitting, 

Fap{z) = ^ + . (2.24) 



Indeed, the derivation given above is one way to derive the Altarelli-Parisi function. We 
illustrate how F{z,y,ai,bi) at 6; = approaches F\p{z) in Fig. 3. 



®Note that there are many different ways to define a momentum fraction variable. The value of the 
splitting function for a given choice of daughter parton momenta does not depend on the momentum 
fraction variable that one uses to label these momenta. We have taken a simple definition of z in order to 
display results in a graph. 
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0.0 0.2 0.4 0.6 0.8 1.0 



z 

Figure 3: The spin averaged splitting function F defined in Eq. (2.23) for a final state g ^ (7 + g 
splitting, plotted versus the momentum fraction z of the gluon, as defined in Eq. (2.20). The 
quark is taken to be massless and we set a; = 4. The four curves are for, from bottom to top, 
y = 0.03,0.01,0.001,0. For ?/ = 0, the result is the Altarelli-Parisi function, [1 + (1 - zf\lz. 

2.2 Initial state q ^ q + g splitting 

Here we consider an initial state q ^ q + g splitting, as illustrated in Fig. 4. For notational 
convenience, we let it be parton "a" that splits, so / = a. We allow both parton "a" and 
parton "b" to have masses, ?na = ?7i(/a) and ?Tib = m(/b). One could, of course, choose 
these masses to be zero. Parton m+1 is a (massless) gluon. The shower evolution for initial 
state particles runs backwards in physical time. Parton "a", which carries momentum pa 
into the hard interaction, splits into the final state gluon with momentum Pm+i and an 
initial state parton that carries momentum into the splitting. For a nearly collinear 
splitting. Pa ~ Va. — Pm+i- In physical time, it is the initial state parton with momentum 
Pa that splits. 




Figure 4: An initial state q + g splitting. 
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Following Ref. [1], we define the kinematics using lightlike vectors pA and pB that are 
lightlike approximations to the momenta of hadrons A and B, respectively, with 2pA'PB = s. 
The momenta of the partons that enter the hard scattering, and p\„ are defined using 
momentum fractions r/a and rj^. After the splitting, the momentum fractions are ?7a and r/b- 
Because parton "a" splits, fji^ ^ rji^. However, with our kinematics, the momentum fraction 
of parton "b" remains unchanged: ?7b = The initial state parton momenta are defined 
to be 

2 

TTT-a 

Pa. = VaPA H PB , 

Ph = 'nhPB-\ -PA , (2.25) 

2 

Tn 

Pa. = fjaPA + ^PB ■ 
VaS 

The momentum of the final state spectator partons changes in order to make some momen- 
tum available to allow both pa and pa to be on shell with zero transverse momenta. We 
denote the total momentum of the final state partons before the splitting hy Q = + Ph 
and after the splitting by Q = Pa + Ph- In the splitting function, we make use of a lightlike 
vector TT-a in the plane of pa and Q. With a convenient choice of normalization, = Pb- 
In the following formulas, it will be convenient to define Pa = Pa — Pm+i- 
Using the definition Eq. (2.8) with the splitting amplitudes v^, from Table 1 of Ref. [1], 
we write the spin averaged splitting function as 

47rQ^g 1 ^ 

2(Pa-PB)M2Pa-Pm+l)^ (2.26) 



X -Tr 
4 



Here Dp,y is given in Eq. (2.11). 

The spin averaged splitting function can be simplified. Let us adopt the notation 

$a = ^ , $6 = ^ . (2.27) 

Then the result can conveniently be displayed in terms of the dimensionless function 

F h_h}±l TFaa . (2.28) 

47ras 

The result is 

p p , Pm+l • PB . ^a^b (??a " Vaf ^ 

^ — -Tcikonal H 1 7Z ^ , ^ — 5- -C^eikonal • (Z.Zy) 

Pa-PB (1 - ^>a^'b) r/l 

Here the first term is the simple eikonal approximation for soft gluon emission, 

Pa- D{pm+1,Q) ■ Pa „„^ 

eikonal — ; ; • [/.6U) 

Pa ■ Pm+l 
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The second term is present in the case of massless or massive quarks and is manifestly 
positive. The third term is present only if and mb are both non-zero. It is likewise 
manifestly positive as long as $a^b < 1- Thus, as for the final state splitting analyzed in 
the previous section, 

F - Feikonal > . (2.31) 




Figure 5: The spin averaged splitting function [r]s,/f}s)F, with F defined in Eq. (2.28), for an 
initial state q ^ q + g splitting, plotted versus the momentum fraction z of the gluon, as defined 
in Eq. (2.35). All partons are taken to be massless. The four curves are for, from bottom to top, 
y = 0.03,0.01,0.001,0. For y = Q, the result is the AltarcUi-Parisi function, [1 + (1 - zY]/z. 

Let us now specialize to ma = mb = and examine the behavior of F in more detail. 
We define a virtuality variable 

Prn+l ' Pa /„ qoN 

y = — (2.32) 

Pa.-Ph 

and a variable representing the momentum fraction of the daughter gluon 

z = P-^^-}P^-Pm+i) (2.33) 

Pa • (Pb -Pm+l) 

We can write z in a different form by using the kinematic relation that is used to define 
the momentum mapping T^a, (Pa + Ph)"^ = (Pa + Pb — Pm+i)^, which is equivalent to 

Pa. ■ Pm+l = {Pa - Pa) ■ Ph ■ (2.34) 

This relation gives 

z = - ^ , 2.35 

1-2/ 



where 



X = 1 - ^ . (2.36) 

Va 
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Note that z and x are equivalent when y = but z varies in the range < z < 1. The 
inverse relation is 

x = z^y{\-z) . (2.37) 



A simple calculation gives 



F = Feikonal + rr^ " V , (2.38) 
1 — X 



where ^ ^ 

-^eikonal — ~ o • (2.39) 
X 

As expected, (1 — x)F = {rig,/f}g)F approaches the Altarelli-Parisi splitting function, [1 + 
(1 — z)"^]/ z as y ^ 0. The approach to the limit is depicted in Fig. 5. 

2.3 Final state g — > g + g splitting 




Figure 6: A final state g ^ g + g splitting. 

Next we consider a final state g — > gg splitting, as illustrated in Fig. 6. According 
Ref. [1], the splitting amplitude is built from the ggg vertex, 

v'^'^'^{pa,Pb,Pc) =v'^'^'^{pa,Pb,Pc) + V2'^^{pa,Pb,Pc)+V^'^'^{pa,Pb,Pc) , (2.40) 

where 

v"^'^{pa,Pb,Pc) = 9"'^{Pa-Pby , 

vf^{Pa,Pb,Pc) = /^(P6 - Per , (2-41) 
V^'^'^{Pa,Pb,Pc) = g^^iPc - Pa^ ■ 

We use v°'^'^ to define the splitting amplitude 

Vl{{p, f}m+l,Sm+l, Sl, Sl) 



^^^"'^ £a{pm+i,Sm+i;Q)*ei3{pi,si;Q)*e''{pi,si;Q) (2.42) 



2pm+l -pi 

X V°''^'^{pm+l,PU -Pm+1 - Pi) D^uiPl + Pm+l] Ui) 
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We have the ggg vertex, polarization vectors for the external particles, and a propagator 
D / {2pm+i-pi) for the ofF-shell gluon. The numerator D^^{pi + pm+i', ?^z) projects on to the 
physical polarization states for the off-shell gluon, 

D^,{P,ni) = -g^, + kr . (2.43) 

P ■ m 



Here n/ is a lightlike vector in the plane of pi and Q, defined as in Eq. (2.5). Then 
njD^i, = 0. 

Following Ref. [1], we define the spin averaged splitting function using Eq. (2.8), 



\m+l-, Sm+l, Si, Si) 



Sl,Sm + l,Sl . 



2 

(2.44) 



+ \v2,l{{P, f}m.+l, Sm.+i, Sl, Sl)\^ - \v3^l{{p, f}m+l, Sm+1, Si, Sl)\" 

Here V2,i and ^3^; are defined as in Eq. (2.42), but with V2^'^ or v^^"^ , respectively, replacing 
the full ggg vertex t>"^'^. Note first of all the prefactor 1/2, which is a statistical factor 
for having two identical final state particles in a g — > g + g splitting. This is the factor Si 
in Eq. (2.8). Then we add |f2,zp — |f3^/p. This does not change the result when we add 
this function to the same function with the roles of the two daughter gluons interchanged. 
With this modification, there is a singularity when daughter gluon m + 1 becomes soft but 
not when daughter gluon I becomes soft. 

One can evaluate Wu as given in Eq. (2.44) by using 

Y,ei.{k,s;Qye,{k,s;Q) = D^,{k,Q) , (2.45) 

s 

where D^y{k,Q) is defined in Eq. (2.11). One might expect a complicated result, but Wu 
is actually quite simple. As in previous subsections, we display the result in terms of the 
dimensionless function 

F = PlA^ Wu . (2.46) 

The result is 

(-A;2)[l + (1- A)2l 

F = i^eikonal + ^ . — ■ 2.47 

4pi-pm+i 

Here i^cikonai is the standard eikonal function given in Eq. (2.14) and 

k^^ = D{pi,Qr,p';,^, , 

PvQpm+i-Q 

Since k_\_, the part of pm+i orthogonal to pi and Q, is spacelike, we again find 

F - Feikonal > . (2.49) 
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Figure 7: The spin averaged splitting function F defined in Eq. (2.46) for a final state g — > g + g 
splitting, plotted versus the momentum fraction z of the gluon, as defined in Eq. (2.20). We set 
ai = 4. The four curves are for. from bottom to top, y ~ 0.03,0.01,0.001,0. For y = 0, the result 
is 2(1 — z)/z + z(l — z). The sum of this and the same function with z ^ 1 — z is the standard 
Altarelli-Parisi function, 2(1 - z) / z + 2z/(l - z) + 2z(l - z). 



We can evaluate F as a function of the variables y and z and the parameter a; , defined 
as for a final state quark splitting in Eqs. (2.17) and (2.20). We find 

1 + fl _ A)2 

F = Feikonal + z{l - z) . (2.50) 

Here i^eikonal was given in terms of z and y in Eq. (2.22) and 

^ = ^. ^T7? vJ ' (2-51) 

x(l -x)(l + y)2 ^ ' 

where the auxiliary momentum fraction x was given in terms of z and y in Eq. (2.21). 
For y ^ 0, F becomes 

F^f{z)^'^^^^ + z{l-z) . (2.52) 



The standard Altarelli-Parisi function, 

fAp{z) = 2 



+ — + z[l-z) 



(2.53) 



is f{z) + /(I — z). Recall from Eq. (2.44) that we broke the symmetry in a g — > g + g 
splitting in such a way that there is a leading singularity for gluon m + 1 becomming soft 
but not for gluon I becoming soft. We could have accomplished the same end by using the 
full ggg vertex but multiplying the splitting function by 6{z < 1/2). Had we done that, 
the small y limit of F would have been f(z) = fx-p{z)9{z < 1/2). This would also give 
f{z) + f{l-z) = fAp{z). 

The full function F{z, y, a;) approaches f{z) as y — > 0, as illustrated in Fig. 7. 
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2.4 Initial state g — > g + g splitting 




Figure 8: An initial state g ^ g + g splitting. 

We now consider an initial state g — > gg splitting, as illustrated in Fig. 8. According 
Ref. [1], the splitting amplitude is again built from the ggg vertex, v'^^'^, Eq. (2.40). We use 
iq define the splitting amplitude for the splitting of one of the initial state partons, 
say parton "a," 



Vi,{{p, f}m+l,Sm+l, Sl, Si) 



2pm+l-Ps. 

X ■u"^'^(pm,+l, -Pa,Pa - Pm+l) D.yu{pa - Pm+V, ria) . 



(2.54) 



We have the ggg vertex, polarization vectors for the external particles, and a propagator 
D / {2pm+i-Ps) for the off-shell gluon. The numerator Z)^;y(pa — Pm+i; ^^^a) projects on to the 
physical polarization states for the off-shell gluon. It is defined using Eq. (2.43), with the 
lightlike vector n.^ = Pb- Following Ref. [1], we use Eq. (2.8) to define the spin averaged 
splitting function from the square of fa. 



i({/;P}m+l) — ^ \va.{{P, f}m+l, ^m+li Sa; '5aJ 



(2.55) 



a, j-^m + l j-^a 



Remarkably, VFaa is rather simple. As in previous subsections, we display the result 
in terms of the dimensionless function 



47ra= 



(2.56) 



The result is 



F — -^eikonal ~l~ 



Pa-Pm+1 



Pm+l-ria 



(Pa -Pni+l)-ns 



+ 



Pa-^a Pm+l-Q +Pm+l-n-a Pa.-Q 
{Pa. - Pm+l ) ■ Pm+l ' Q 



(2.57) 
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Here -Feikonai is the eikonal function, Eq. (2.30), and /c^ = D{pi^, Q)^u Pm+i ^Iq. (2.48). 
Examination of Eq. (2.57) shows that, as in the previous cases. 



eikonal 



> 



(2.58) 



To see this, one needs to know that spUtting kinematics ensures that (pa — Pm+i)"f^a > 0. 
We note that the sphtting kinematics allows non-zero parton masses, although the gluon 
that splits here is massless. 

Let us look at this assuming massless partons and using the splitting variables y, z 
and X = z + y{l — z) defined in Sec. (2.2). A straightforward calculation gives 



X X 



% + 2{l -y)z 



(l-?/)2 



+ 



2z{l-y)+y 



l-{l-y)z) z{l-z){l-yY + y 



(2.59) 



As expected, {l—x)F = {r]s,/fi^F approaches the Altarelli-Parisi splitting function, 22;/(l- 
z) + 2(1 — z)/ z + 2z(l — z) as y — > 0. The approach to the limit is depicted in Fig. 9. 



200 




Figure 9: The spin averaged splitting function {i]^/fi^)F, with F defined in Eq. (2.56), for an 
initial state g ^ g + g splitting, plotted versus the momentum fraction z of the gluon, as defined 
in Eq. (2.20). The four curves are for, from bottom to top, y ~ 0.03,0.01,0.001,0. For y — 0, the 
result is the standard Altarelli-Parisi function, 2(1 — z)/z + 2z/{l — z) + 2z{l — z). 



2.5 Other cases 

We have covered the cases of quark or gluon splittings in which a daughter gluon enters the 
final state. There is also the possibility of an antiquark splitting replacing a quark spitting, 
but, because of charge conjugation invariance, these are essentially the same as the quark 
splitting cases. There are also cases in which no daughter gluon enters the final state: final 
state and initial state g — > Q + Q and initial state q ^ q + g and g — > g + g in which the 
gluon enters the hard scattering and the quark or antiquark enters the final state. The spin 
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averaged splitting functions for these cases are manifestly positive. In these cases, there is 
no leading singularity when a final state daughter parton becomes soft, so we do not need 
to consider soft gluon singularities. We list the results for these cases in Appendix A. 

3. Interference diagrams 




Figure 10: An interference diagram. A gluon, labelled m + 1, is emitted from parton I in the 
amplitude and from parton k in the complex conjugate amplitude. This diagram has a leading 
singularity when the gluon is soft. 

We have analyzed the spin averaged splitting functions Wu, which correspond to the 
squared amplitude for a parton I to split into daughter partons with labels / and m+1. Now 
we need to consider interference diagrams, such as the diagram illustrated in Fig. 10. In 
the amplitude, parton / can change into a daughter parton with label / by emitting a gluon 
with label m + 1. In the complex conjugate amplitude, parton k can change into a daughter 
parton with label k by emitting a gluon with label m+1. If we were to temporarily ignore 
questions about how to define the kinematics and were to use the splitting amplitudes vi 
and Vk for this, the corresponding contribution to the splitting function would be 

Vl{{p, f}m+l,Sm+l,Sl,Sl)5s,^,s,, Vk{{p, f}m+l, s'm+l, s'k, s'k)* 5^^s'^ . (3.1) 

This function is singular when gluon m + 1 is soft, Pm+i — > 0. However it does not have 
a leading singularity when gluon m + 1 is collinear with parton / or parton k. For this 
reason, we can use a simple eikonal approximation to the splitting amplitude, 

soft/r- r-i ^ X r-, e £{Pm.+l, Sm+l] Q)* 'Pi n\ 

Vl [{P, f}m+l,Sm+l, Si, Si) = V^T^as Ss^^si 1 ^ , (3.2) 

Pm+I-Pl 

if parton m + 1 is a gluon, with vf°^^ = otherwise. 



- 16 - 



Making the eikonal approximation, the sphtting function is 



(3.3) 



This function gives the dependence of the sphtting operator on momentum and spin. In 
Ref. [1] there is a separate factor that gives the color dependence. This factor is an operator 
on the color space that we can call tj i8>t^, where t| is the operator in Eq. (2.2) that inserts 
the proper color matrix into line / in the amplitude and t/^ inserts the proper color matrix 
into line k in the complex conjugate amplitude.^ We do not yet make any approximations 
with respect to color. In Sec. 5, we will make the approximation of keeping only the leading 
color conributions. 

There is an ambiguity with the prescription (3.3). The functions in Wik are defined 
in terms of daughter parton momenta and flavors, {p, f}m+i- However, we want to define 
{p, f}m+i from the momenta and flavors {p, f}m before the splitting together with the 
splitting variables {CpjCf}- We need to specify what relation to use. One way is to use 
the kinematic functions that we use for the splitting of parton I into a daughter parton I 
and the gluon m + 1, {p, f}m+i = Ri{{p, f}m,{CpXi})- With this mapping, we define a 
function wj^j^ of the {p, f}m and the splitting variables. Alternatively, we could use the 
kinematic functions, Rk, that we use for the splitting of parton k into a daughter parton 

(k) 

k and the gluon m + 1. With this momentum mapping, we define a function W^j, of the 
{Pj f}m and the splitting variables. Instead of using one or the other of these possibilities, 
we average over them. We use wj^l^ with weight Aik and wj^^^ with weight A^i- In Ref. [1], 
we let the weight functions take the default value Aik = A^i = 1/2. This choice is certainly 
conceptually simple. However, we can obtain spin-summed splitting functions that have 
nicer properties if we define the weights as certain functions ^;fc({p}m+i) and ylfc;({p}m+i) 
of the momenta. It is simplest to specify the functional forms of the weight functions using 
the momenta {p}m+i after splitting. The momenta after splitting are to be determined 
by the mapping TZi for An^ and by the mapping 7^^, for Af^i.^ The weight functions are 
non-negative and obey ^/fc({p}m+i) + Aki{{p}m+i) = 1 at fixed momenta {p}m+i- The 
relation Aik+A^i ~ 1 then holds at fixed {p, f}m and splitting variables. This approximate 
relation becomes exact in the limit of an infinitely soft splitting, for which the mappings 
TZi and TZk become identical. 

With the choice of momentum mappings determined by Aik and A^i, the net splitting 
function, including the color factor, summed over the two graphs arising from interference 
of soft gluons emitted from partons / and k, is 



AikW^^ + AkiWi'^ 



4 + 



A 



Ik 



Ik ''I 



AkiWtf^+Ai,wtf 



ti 



(fc) 4 



ik + w;:r tl ® t 



(3.4) 



^In Ref. [1], we write tj{fi fi +g) for the operator that we here call just tj and we denote the operator 
tj(g)tf, hy g{i,k). 

*This is expressed most precisely using the operator language of Eq. (8.26) of Ref. [1]. 
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We will see in the following section that we obtain spin-summed splitting functions 
that have nice properties if we define Ai^ as a ratio, 



Aiki{p} 



m+l) 



Bikm 



m+l J 



Blki{p}m+l) +Bkli{p} 



m+l) 



where 



Bik{{p}m+i) = ^I^^^ prD{pm+i,Q)-pi . 



Pm+l -pi 

Here D{pm+i,Q) is the transverse projection tensor defined in Eq. (2.11). 



(3.5) 
(3.6) 



4. Spin-averaged interference graph splitting functions 

The part of the soft splitting function representing l-k interference that is associated with 
the kinematic mapping TZi is 



A, 



Ik 



(4.1) 



We now make the approximation of setting s' = s for each parton, summing over the 
daughter parton spins, and averaging over the mother parton spins. The sum over spins of 
W/P is the same as the sum over spins of TV^P. Thus the spin averaged splitting function, 
including the color factor, becomes^ 

1 
2 



tk+i 



t 



{s'}={s} 



(4.2) 
(4.3) 



where 

Wik = ^ Aik 
Here we have used the notation {s'} = {s} to indicate the instruction to set s[ = si, s'l = si, 

k 



s'l. = Sfc, s'f^ = Ski and s'j^^j^i = Sm+i- The structure of Wik is quite simple 



Wik = 47ras 



Ik 



Pi ■ D{pm+i;Q) -Pk 

Pm+l- Pi Pm+l-Pk 



(4.4) 



We can associate Wik with the splitting of parton I, since it uses the kinematic mapping 
TZi- Then we are led to consider the relation of Wi^ to the direct splitting function Wu. 
Now, the color factor that multiplies Wu is tj (E> t^. However, as discussed in Ref. [1], the 
invariance of the matrix element under color rotations implies that 



kf^l 



h + i 



(4.5) 



Thus we can combine the direct and interference graphs to give 



*fc + ^k 



[Wu - Wik] 



(4.6) 



^The function Wik here equals the product 2An:Wik of functions in Ref. [1]. 
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We will see in Sec. 5 that the color factor here is very simple in the leading color limit, 
essentially amounting to multiplying by Cp or zero. We are thus motivated to investigate 
the coefficient of this color operator, Wu — Wik- 
It is useful to break Wu — Wik into two pieces, 

Wn - Wiu = {Wu - °°'^') + (Wf ""'^^ - Wiu) ■ (4.7) 
Here, we recall from Eq. (2.10), 



Wu = 47ras • 



We have investigated {Wu — Wu ) in Sec. 2 and found that 

Wu - Wf" > . (4.9) 
Thus we should consider — Wik- We have 

p^eikonal _ ^ ^ A-Ka^ | Pi ■ D{pm+i;Q) 'Pi _ ^^^^ Pi ' D{pm+l]Q) ' Pk 1 ^^-^ 
" Pm+I-Pl \ Pni+l-Pl Pm+l-Pk J 

We can simplify this if we use the definitions (3.5) of Ai^ and (3.6) of 



-ppy cikonal ^^yt 

Wu - Wik 
47ras 



1 ipm+l -Pk . . „ , . n " \ 

— \ Pi- D -pi- 2Aik pi- D -pk) 

i-Pki Pm+i ■ Pi J 



Pm+l ■ Pi Pm+1 ■ Pk I Pm+1 ' Pi 

47ras 1 



Pm+1 • Pi Pm+l ■ Pk 

47r«s Afc 
Pm+l ■ Pi Pm+l ■ Pk 



Bik - '^Aik Pr D-pk 



Bik + Bki-2pi-D-pk 



AtTUs Alk f Pm+l -Pk ^ n - , Pm+l ' Pi . ^ - 

Pi- D ■pi + Pk- D -pk 



Pm+l - Pi Pm+l - Pk [ Pm+l - Pi Pm+l - Pk (^-H) 

-2 pi- D -pk 

= — — — {{Pm+l - Pk) Pi- D - pi 

Pm+l - Pi Pm+l - Pk Pm+l - Pi Pm+l - Pk 

+ {Pm+l - Plf Pk- D-pk-2 {pm+l - Pk){Pm+l - Pl) Pl - D - Pk] 

_ Alias Alk 1 

Pm+l - Pl Pm+l - Pk Pm+l - Pl Pm+l - Pk 

X {Pm+l - Pk Pl- Pm+l - Pl Pk)- D - {pm+l - Pk Pl - Pm+l - Pl Pk) ■ 

We can simplify this further by noting that the vector Pm+i -Pk Pi —Pm+i -Pi Pk is orthogonal 
to Pm+l, so that only the term —g^'^ in D^'^ contributes. Thus 

^eikonal ^ AtTOs Ak -{Pm+l - Pk Pl - Pm+l - Pl Pk)'^ 

^11 - ^Ik = -z- ; ; 1. • (,4.i2j 

Pm+l - Pl Pm+l - Pk Pm+l - Pl Pm+l - Pk 
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Since the vector Pm+i ■ Pk Pi — Pm+i • Pi Pk is orthogonal to the Ughthke vector Pm+i, it is 
either Hghthke or spacehke. Furthermore, Aik > 0. Thus 

^eikonal _ > q _ (4_13) 

Thus both parts of our sphtting function, Wu — 1^;"^^°°^^ and Wil^°^^^ — Wik, are non- 
negative. This means that we can use these functions as probabihties in constructing a 
parton shower Monte Carlo program without needing separate weight functions. We discuss 
this further in Sec. 6. 

The analysis so far has allowed partons / and k to have non-zero masses. Let us now 
consider the case of massless partons, pf = p\ = The massless result can be understood 
in more detail if we write it in terms of three- vectors in the frame in which Q = 0. We 
define Uk-, ui, and Um+i to be unit three- vectors in the directions of the space parts of pi, 
Pk, and Pm+i respectively. Then 

TTjeikonal TTr 4:TTO:s 2Q^ / 1 -, i\ 

Wu -Wik = -^^ — — — g[u„,+i,ui,Uk) , (4.14) 

[Q-Pni+l)^ (1 - Um+l-Ul) 

where 

g{um+i,ui,Uk) = J- = -.w-i I ^ ^TTTi = \^^ , ^ (4-15) 

(1 - Um+1-Ul){l + Um.+l-Uk) + (1 - Um+rUfc)(l + Um+l-Ul) 

We can make some comments about this. First, the splitting probability is singular when 
the angle between Um+i and ui approaches zero, (1 — Um+i ■ ui) — > 0. This is the standard 
coUinear singularity, seen in the soft limit. Second, when (1 — Um+i -ui) <C {1 — Uk -ui) <^ 1, 
Wu — Wik behaves like 1/(1 — Um+i • ui). If we integrate over the angle of Um+i with 
a lower cutoff on the angle between Um+i and ui, the integral is logarithmically sensitive 
to the cutoff. Third, when (1 — -u^ • ui) ^ (1 — Um+i • ui) <^ 1, Wu — Wik behaves 
like 1/(1 — Um+i • ui)"^. If we were to put an upper cutoff on the angular integration, there 
would be no logarithmic sensitivity to this cutoff. Thus, only the angle ordered region 
(1 — Um+i ■ ui) ^ (1 — "Ufc • ui) is important in the integral over angles. There is a smooth 
decrease in the splitting probability when the angle between Um+i and ui becomes greater 
than the angle between Uk and ui. There is no sharp cutoff. 

We illustrate this in Fig. 11. We take the polar angles of Um.+i to be 6,(j) where ui 
is along the 9 = axis. We choose Uk to have polar angles 9k = 0.1 and (j)k = 0. Then 
we plot g{um+i,'Ui,Uk) versus 9x = 9cos(j) and 9y = 9 sin (p. Since Wu — Wik oc 
g{um+i,ui,Uk)/i^ — Um+i • ui), the main feature of Wu — Wik is a singularity at 
9x = 9y = 0. We see that the factor g that multiplies the singular factor is a smooth function 
with a gentle peak between ui and Uk- This peak above g = 1 represents constructive 
interference. When Um+i moves outside the "angle ordered" region (1 — Um+i • ui) < 
(1 — Um+i • Uk), the factor g drops below 1 and decreases to zero, representing destructive 
interference. We notice in Fig. 11 that there is an enhancement of soft gluon radiation in 
the region between the directions of parton I and parton k. This enhancement is known as 
the string effect and has been observed experimentally [7]. 
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0.1 



Figure 11: The function g defined in Eq. (4.f 5) that serves to suppress soft gluon radiation outside 
of the "angle ordered" region. The plot coordinates are 9x = Ocoscfi and 9y = dsm(p, where 0,0 
are the polar angles of Um+i- The vector is at = and the vector Uk is at 9 ^ OA, (f) ~ 0. 



5. The leading color limit 

We have studied the spin-averaged splitting function [Wn — Wik]- Here Wn describes the 
square of the graph for emission of a gluon from parton /. There are also interference graphs 
between emitting the gluon from parton I and emitting the same gluon from parton k. The 
function Wi^ describes the part of the interference graphs that we group with parton /. 
These functions give the momentum dependence. They multiply a color operator as given 
in Eq. (4.6), 

(5.1) 



1 r 
2 



+ A ® h 



We have so far not made any approximations with respect to color. Let us now take 
the leading color approximation. To do that, recall from Ref. [1] that we use color states 
based on color string configurations. For instance, we could have a state [4, 5, 2, 3, 1] in 
which 4 labels a quark, 1 labels an antiquark, and 5, 2, and 3 label gluons that connect, in 
that order, to a color string between the quark and antiquark. One can also have a closed 
string such as (4,5,2,3, 1) in which all of the partons are gluons. A color basis state can 
also consist of more than one color string connecting the partons. In general, the amplitude 
can have one color state |c) and the complex conjugate amplitude can have a color state 
^c'l with c' 7^ c. However, in the leading color approximation we can only have c' = c. 
Additionally, in the leading color approximation we have 

-t] ® = -4 = -i [4 ^tk + tl(S) tij ~ Cf al a,fc . (5.2) 

Here represents the operator that inserts gluon m + 1 between partons / and k on the 
color string if these partons are adjacent to each other on the same color string, that is, if 
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partons / and k are color connected. When aj^ is applied to a state |c) in which / and k are 
not color connected, we define a|^|c) = O.^'^ For the complex conjugate amplitude, (c| aik 
again gives a state with the soft gluon inserted between partons / and k. Thus, starting 
with a color state |c) in the amplitude and (c| in the complex conjugate amplitude, we 
get zero if partons / and k are not color connected and we get a new color state with the 
soft gluon inserted between / and k if / and k are color connected. The bookkeeping on 
color connections is a standard part of parton shower event generators. The momentum 
dependent numerical factor [Wu — Wik] is multiplied by a color factor Cp. 

This analysis has covered the case in which parton m + 1 is a gluon, so that there are 
interference graphs arising from this gluon being emitted from parton / in the amplitude 
and from parton k in the complex conjugate amplitude (or the other way around). There 
are also graphs for which parton m + 1 is a quark or antiquark, as described in Sec. 2.5. 
In these cases, we have just the splitting function Wu, which multiplies the color operator 
t\ ^t^. This operator is very simple in the leading color limit. 

Consider first the case of an initial state splitting in which fi = q and {fi, fm+i} = 
{g, q}, where g is a quark flavor (u, u, d, . . . ). In physical time, this is a splitting g — > 5 + 
while in shower time it is a splitting g — > g + g. As discussed in Sec. 7.3 of Ref. [1], 

tl(E)ti = Cf (E) ag{l) . (5.3) 

Here ag{l) represents the operator that inserts the gluon at the end of the string terminated 
by quark I before the splitting and terminated by quark m+1 after the splitting. -"^^ Similarly, 
in the case of an initial state splitting in which fi = q and {fi,fm+i} = {gj^ji we have 
the same result, where now a\{l) represents the operator that inserts the gluon at the end 
of the string terminated by antiquark I before the splitting and terminated by antiquark 
m + 1 after the splitting. 

Consider next the case of an initial state splitting in which fi=g and {fi,fm+i} = 
{q, q}. In physical time, this is a splitting q ^ q + g, while in shower time it is a splitting 
g ^ q + q. As discussed in Sec. 7.3 of Ref. [1], in the leading color limit, 

4 ® ~ Tr 4(0 ® ag(/) , (5.4) 

where Tr = 1/2 and aq{l) splits the color string at the point at which gluon / attaches, 
creating new string ends corresponding to the quark and the antiquark. The same analysis 
applies for an initial state splitting with fi = g and {fi, fm+i} = {q, q} and for a final state 
g ^ g + ^ splitting, for which {fi, fm+i} = {q, q}- 

6. Evolution equation 

We now have the information that we need to present the formulas from Ref. [1] for parton 
shower evolution specialized to the spin averaged, leading color approximation. In the 

^°We here adapt the notation of Ref. [1], where we had gluon insertion operators aJj_(Z) and a^_{l) that 
insert the gluon to the right or the left of parton I, respectively. If partons / and k are color connected, 
we have a]^j.|c) — a^(Z)|c) or alf,\c) — a'_{l)\cj, depending on whether parton k was to the right or left of 
parton I along the string. 

^^This operator is denoted a^(i) in Ref. [1]. 



- 22 - 



general case, we had basis states \{p, f,s',c',s,c}m) with two color configurations {c}m 
and {c'}m, representing the color state in the amplitude and the color state in the complex 
conjugate amplitude, respectively, and two spin color configurations {s}m and {s'}m- In 
this paper, we have averaged over spins, so that we can describe the evolution of the states 
without referring to spin at all. We also use the leading color approximation, so that we 
always work with states with {c}m = {c'}m- Thus our description is vastly simplified and 
we can work with basis states \{p, /, c}m)- 

As in Ref. [1], we use the logarithm of the virtuality of a splitting as the evolution 
variable, so that a splitting of parton I is assigned to a shower time t = Ti{{p, 

mpju+i) = iog( — - — — ) , (6.1) 

where fi = fi+fm+i and Qq is the starting virtuality scale. Shower evolution is based on the 
probability that, at shower time t, a state \{p, f, c}m) that had not already split now splits 
to make a new state \{p, f,c}m+i) with one more parton. This probability is represented 
as a matrix element of a splitting operator Hf'\t), which is similar to the splitting operator 
7ii{t) of Ref. [1] except that the spin averaged, leading color approximations ("(0)") have 
been applied. Then Tl[^\t) operates on states \{p, f, c}m) instead of the states of the full 
theory. We write 

{{p,f,cU+i\n[^\t)\{p,f,cU) 

E nc{a)nc{b)i]a,rib fa/A{f]!^,lJ^F)fb/B^'^^ 
/ j-f^ _L J _ 1 

^ nc(a)nc(fe) fj^fjy, fa/AiV^.-, I^F)fb/B{r]h,fJ'F) 

rn+l 



X !.9{frn+l=g) ^ {{c}rn+l\alk\i^}rn)'^lk{{P,f}rn+l) ^ 
it 

+ 9{fm+l^g)e{fl=g) {{cU+i\al{l)\{cU) '^U{{pj}m.+l) 

+ 9{L+i^g)e{fi=g) {{cy^+i\al{l)\{cU) '^u{{pJU+i] 

The first line on the right hand side of this formula contains factors copied directly from 
Ref. [1]. There is a sum over the index / of the parton that splits. Then there is a ratio 
of parton distribution functions. This ratio is 1 for a final state splitting but different 
from 1 for an initial state splitting. The next line concerns the relation of the variables 
{p, /}m+i and t to the variables {p, f}m- For the flavors, this factor vanishes unless there 
is a QCD vertex for fi^fi + fm+i and it vanishes unless fj = fj for the other partons. 
For an allowed relationship between {/}m+i and {f}mi the flavor factor is 1. There is 
a similar factor for the momenta. Given the momenta {p}m, the momenta {p}m+i must 
lie on a certain three dimensional surface specified by the momentum mapping TZi defined 
in Ref. [1]. The function {{p, f}m+i\'Pi\{p, f}m) contains a delta function on this sur- 
face. There is also a delta function that defines the shower time t. Thus if we integrate 
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({p, /, c}m+i|Wj /, c}m) over t and the momenta {p}m+i, we are really integrating 

over three variables that describe the splitting of parton I. 

The final factor in Eq. (6.2) contains three terms. Our main interest is in the first term, 
for fm+i = g- There is a sum over the index k of other partons in the process. These are the 
partons that might be connected with parton I in an interference diagram. The remaining 
factors are rather complicated in the general case described in Ref. [1] , but are quite simple 
in the spin averaged, leading color approximation. The factor ^{c}m+i|a[fc|{c}m) embodies 
the color considerations described in Sec. 5. It equals 1 provided two conditions hold. First, 
partons / and k must be color connected in the initial color state {c}m- Second, the new 
color state {c}m+i must be the same as {c}m with the gluon with label m + 1 inserted 
between partons / and k. If either of these conditions fails, this factor vanishes. The 
remaining factor is the splitting function 

^^^ = CF[Wa-Wik] . (6.3) 

We have seen explicitly what this factor is, and have noted that is positive. 

The next term in the braces in Eq. (6.2) applies to an initial state splitting in which 
fl = g and {fi,frn+i} is either {q,q} or {q,q}. The color factor ({c}m+i|4(/)|{c}m) is 
1 if the new color state {c}m+i is the same as {c}m with the end of the string at quark 
or antiquark / now terminated at quark or antiquark m + 1 and the new the gluon with 
label I inserted just next to the end of the string. Otherwise, this factor vanishes. The 
corresponding splitting function is 

$zz = CfWz, . (6.4) 

The final term in the braces in Eq. (6.2) applies to an initial state splitting in which fi = g 
and {fi,fm+i} is either {q,q} or {q,q}. The color factor <|{c}m+i|aq(0|{c}m) is 1 if the 
color state {c}m+i is related to {c}m by cutting the color string on which parton / (a gluon) 
lies into two strings, terminating at the new quark and antiquark. Otherwise, this factor 
vanishes. The corresponding splitting function is 

^ii=TnWii . (6.5) 

We have now specified the probability that a state \{p, f, c}m) splits. The probability 
that this state does not split between shower times t and t' is 

A(0)(t,t';{p,/,c}„)=exp(^-^^*dT (l|wf (T)|fe/,c}^)^ . (6.6) 

Here [l\Tlf'\T)\{p, f, c}m) is the inclusive probability for the state \{p, f, c}m) to split at 
time r, 

(l|<)(r)|{p,/,c}^) = I [d{pJ,cU+i] {{pJ,cU+i\n['\t)\{p,f,cU) . 

(6.7) 
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To get the inclusive splitting probability, we have integrated over the momenta {p}m+i 
after the splitting and summed over the flavors and colors, using the integration measure 
in Eq, (3.15) of Ref. [1], supplemented by a sum over color states. 

With these ingredients, we can describe shower evolution using the evolution equation 
(14.1) from Ref. [1]. The evolution from a shower time t' to a final time tf at which 
showering is terminated is given by an operator U^^\t{,t') that obeys^^ 

U^''\ti,t')=M^'\tut')+ j\TU^'\tuT)n['\T)M^'\T,t') . (6.8) 

Here Af^'^\t',t) is a no-splitting operator defined by 

AA(0)(t',i)|{p,/,c}^) =AW(t',t;{p,/,c}^)|{p,/,c}„) . (6.9) 

If we apply this to a state \{p, /, c}m) that exists at shower time t', we have 

(tf , t') I {p, f, cU) = A(0) (tf , t'; {p, f, cU) I {p, f, c}m) 

+ ^/^r ^^1^^, I [d{pJ,cU+i] Z^(0)(tf,r)|{p,/,cWi) 

X ({p,/,cWi|<)(r)|{p,/,c}^)A(0)(r,t';{p,/,cW . 

(6.10) 

The first term gives the probability that the state does not split before shower time tf . The 
main evolution is represented by the second term. There is an integration over the shower 
time r of the next splitting and over the splitting parameters. In an implementation of 
this equation, the integration would be performed by Monte Carlo integration. That is, 
we would choose r and {p, /, c}m+i with some probability density p that contains delta 
functions that restrict r and {p, /,c}m+i to the allowed surface defined by the Eq. (6.1) 
for r and the momentum mapping TZi. Then we multiply by a weight w defined by 

{{pJ,c}m+i\n[''\r)\{p,f,c}m)A^'Hr,t';{pJ,cU) = wxp . (6.11) 

In the present case, the integrand has two welcome features. First, it is positive. Second, 
using the definition of A^''^ 

(^^^Tl)! r^^^ J f"^^^' ^}rn+i\n['\T)\{p, /, c}m) A(°)(r, t'; {p, /, cU) 

= 1 . 

(6.12) 

^^According to Eq. (3.15) of Ref. [1], there is an extra normalization factor ({c}m+i |{c}m+i) in Eq. (6.7). 
With our choice of the normahzation of color states, this factor is not exactly 1, but it is 1 in the leading 
color limit. 

^^In Ref. [1], [TIi{t) — Vs(t)] appears in place of H^^^t) here. With the leading color approximation, 
Vs(r) = 0. 
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Thus the function 

{{p, f, c}m+i I (r) I {p, /, cU) A(0) (r, t'; {p, f, cU) 

P = 4 Tvr^ (6.13) 

^ (m + 1)! ^ ^ 

is positive and properly normalized to be a probability density. Using standard methods 
from shower Monte Carlo algorithms [8, 9, 10, 11] , we can choose points with this probability 
density. Then w = 1. With a probability A^^\ti,t';{p,f,c}m), the point selected will be 
in the range if < r < oo. In this case, there is no splitting and we simply keep the state 
\{p, f,c}m)- This corresponds to the no splitting term in Eq. (6.10). If r < t{, the state 
splits to {p, /,c}m+i. Then, according to Eq. (6.10), we should apply U^^\ti,T) to this 
state, repeating the process. Thus the evolution proceeds by what is known as a Markov 
chain. 

The starting point for evolution is a state that is a mixture of the basis states | {p, /, c}m) 
for m = 2, assuming that we start with a 2 ^ 2 hard process, 

\P^'\^)) = ^ / [d{pJ,c}2] |te/,c}2)({p,/,c}2|/9(°)(0)) . (6.14) 

Here ({p, /, c}2 1/0^'''*(0)) is obtained from the 2 — > 2 matrix element summed over spins, 

To implement Eq. (6.14), one would choose points {p, f, c}2 by Monte Carlo methods. This 
gives the starting point for the shower evolution. The state [/^'^'^^(O)) then evolves into a 
state 

|pW(tf))=Z^W(tf,0)|/5W(0)) (6.16) 

at the shower time tf at which we choose to terminate shower evolution. At this point, 
as described in Ref. [1], the desired cross section is obtained by applying a hadroniza- 
tion model to the component states /, cIat) in \p^^\t{)^, producing a hadronic state 
Z^'^^^(oo,if)|p(0)(tf)). Then the desired cross section (^[-Fh] results from applying the mea- 
surement function to the hadronic states produced. Thus 

^(°^[Fh] = (Fh|^/'^^''(oo,tf)|p(0)(tf)) 

^ / [d{p,f,c}N] {F^\u'^-''{^,ti)\{pJ,c}N){{p,f,c}N\p^''\k)) . 



N 



(6.17) 



Just as in the parton shower evolution, the integration in Eq. (6.17) can be implemented 
by simply taking the states /, cIat) generated by the shower evolution and passing 
them to a Monte Carlo implementation of a hadronization model. Then application of 
the measurement function is acheived by, for instance, putting the events into desired bins 
according to the momenta of the resulting hadrons. 



14 , 



*As explained in Ref. [1], we should most properly project out the component of \M{{p, f}2)) that is 
proportional to a color basis state |{c}2) by using a dual basis state ^({c}2|, but in the leading color limit 
there is no distinction between the dual basis states and the ordinary basis states. 
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7. Other approaches 



In this section, we sketch the relation of the shower evolution of this paper to some other 
approaches to the description of parton showers. For the shake of the simplicity we work 
only with massless partons in this section but it is still allowed for the non-QCD particles 
to have non-zero masses. 

7.1 Dipole shower 

One possibility for organizing the gluon radiation in a (spin averaged, leading color) parton 
shower is to use the same functions that are used for organizing the subtractions in a next- 
to-leading order perturbative calculation. In particular, the dipole subtraction scheme of 
Catani and Seymour [6] is an attractive possibility [12] that has been developed as the 
basis for parton shower programs by Schumann and Krauss [4] and by Dinsdale, Ternick 
and Weinzierl [5]. 

To see how this can work, consider the case that the emitted parton m -|- 1 is a gluon, 
so that the splitting operator is given by the main term in Eq. (6.2), 

{{pJ,c}m+i\n['\t)\{pj,cU) 

I nc(a)?lc(&) ??a?)b fa/AiVs., /^|)/fe/i?(^b, /^|) , . 

X {{P, hm+llVllip, f}m) 6(t - Ti{{p, f}m.+l)^ ({c}m+l | | {c}m) 
X ^lk{{P, f}m+l) ■ 

The term /, k generates gluons predominately soft or collinear with parton I. That is 
because is singular when pm+i is soft or collinear with pi but finite when Pm+i is 
collinear with p^ . Each term is defined with its own phase space mapping Vi and evolution 
parameter t. Now we can use the momentum mappings of Catani and Seymour. These 
obey 

{{pJ}m+l\Vfk\{P,f}m) ~ {{P, f}m+l\Vl\{p, f}m) when p^+i Xpi 

{{p,f}m+i\rf^\{p,fU) ~ {{pJU+il-PilipJU) (7.2) 

~ {{P, hm+ilVklip, f}m} when p^+i ^0 . 

We can also use the splitting functions <I>^| of Catani and Seymour. These substitutions 
give 

{{pJ,cU+i\nr{t)\{p,f,cU) 

_ y^y^ nc(a)nc(fe) r?a?7b /aM(^a, /^|')/6/5(??b, /^l') 

X {{p,f}rn+l\rE\{pJU)s(t-Ti{{pJU+i)'^ {{c}m+l\4k\{^U) 
X ^Tk{{pJ}m.+l) ■ 
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The splitting operator Ti.f^{t) matches TL^\t) in the collinear and soft hmits. 

We see that the structure of shower generation using the Catani-Seymour functions is 
quite similar to that of this paper. It is of interest to compare the splitting functions in 
the soft limit, Pm+i 0. Using the definitions in Ref. [6], we have 

mi{pJ}m+i) ~ 9^^{n^+,,ui,n,;Ei/E,) , (7.4) 

[Q-Pm+l) (1 - Um+l-Ui) 

for Pm+1 0, where 

g^'^Um+i^uu Uk; Ei/Ek) = — ^=rVTTi = ^ ■ C^-^) 

Here Ei and Ef^ are the energies of partons / and k, respectively, in the rest frame of Q, 
the total momentum of the final state partons. Thus Ei/Ek = pi ■ Q/pk • Q- This function 
is similar in form to the function g of this paper, plotted in Fig. 11, but it depends on 
the ratio Ei/Ek- We plot it in Fig. 12 for Ei/Ek = 3 and Ei/Ek = 1/3. We see that the 
Catani-Seymour functions assign little soft radiation to the more energetic of partons I and 
k. More soft radiation is assigned to the less energetic parton of I and k, with quite a lot 
of the radiation going in approximately the direction of the more energetic parton. 

The the final state shower in the latest version (version 8.1) of Pythia [3, 10] is essen- 
tially a dipole shower as described above. In particular, the splitting function describing 
gluon emission in the soft limit Pm+i ^ is that in Eq. (7.4) with the same function g as 
given in Eq. (7.5). 

7.2 Antenna shower 

In the method of this paper and in a dipole shower following the Catani-Seymour scheme, 
the creation of a new gluon is attributed to the splitting of one of the previously existing 
partons. This requires that for the interference graph between the amplitude for emitting 
the gluon from parton I and the amplitude for emitting the gluon from parton k, one assigns 
a certain fraction A of the graph to the splitting of parton I and a fraction 1 — A to the 
splitting of parton k. In an antenna shower, one treats the pair of color connected partons, 
/, A; as a unit. The l,k dipole constitutes an antenna that radiates the daughter gluon. 
The pioneering development along these lines is the final state shower of Ariadne [2]. 
More recent examples include those in Refs. [13, 14]. There is a corresponding subtraction 
scheme for next-to-leading order calculations, antenna subtraction [15]. 

To define an antenna shower, we choose a momentum mapping Vj^^^ with the properties 
previously defined and with the symmetry property 

vr^'=vir . (7.6) 

We also redefine the shower evolution variable to be symmetric under I <-> k interchange. 
For instance, we could take 

t = log(——^^^^^—) . (7.7) 
\2mm[prpm+i, pk-Pm+i\J 



^^One ought to call this a dipole shower, but then one would need a new name for the kind of shower 
described in the previous subsection. 
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Figure 12: The soft radiation function for parton I corresponding to Catani-Seymour dipole 
splitting. The function g defined in Eq. (4.15) and plotted in Fig. 11 is replaced by the function 
^c.s., Eq- (7.5), plotted here. In the plot on the left, Ei/Ek = 3. In the plots on the right, 
Ei/Ek = 1/3. The plot coordinates and value of 6ik are as in Fig. 11. 



Then we can rewrite the sum over / and /c as a sum over pairs /, A;, with each pair counted 
once, giving 

{{pj,c]m+l\nr\t)\{pj,c]m) 



^ (m + 1) 



nc{a)nc{b) f}J}\, fa/Aiv-^, /^f) A/B(??b, fJ-p) 



pairs (' 

Ql 



X {{p, fu+i\rtt\{p, fu) 5[t - log 

X {{c}m+i\al\{cU) (te/}™+i) . 



2min[prPm+i, Pk-Pm+i] 



Here can be 



a>ff = <5ik + <5m (7.9) 



or any function that matches it in the soft and collinear limits. 

In the soft limit, pm+i 0, $f^* approaches the soft limit of the sum + ^ki, which 

is 

There is no function g here. The function g in the previous subsections arises from sepa- 
rating this into two terms, one that remains finite when (1 — Um+i-Uk) — > and the other 
that remains finite when (1 — Um+i-ui) 0. 
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7.3 Angular ordering approximation 
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Figure 13: The angular ordering approximation. The function g defined in Eq. (4.15) and plotted 
in Fig. 11 is replaced by the function g^.o., Eq. (7.11), plotted here. The plot coordinates and value 
of 9ik are as in Fig. 11. 

With massless kinematics, the distribution of soft radiation that is kinematically of 
the form for a splitting of parton I is proportional to g{um+i,ui,Uk)/{l — Um+i • ui), as 
given in Eq. (4.14). From the plot of g in Fig. 11, we see that the soft gluon radiation 
from partons I and k is approximately confined to a cone between pi and pk ■ This is called 
"angular ordering." There is also an angular ordering approximation [16] that is sometimes 
used for parton showers and, in particular, lies at the heart of Herwig [11]. With this 
approximation, the function g in Fig. 11 is approximated by the function plotted in Fig. 13, 



We see that in the angular region between the two hard parton directions {0x ~ 0.5, ~ 
in the figures), the angular distribution of the soft radiation determined by the exact 
function g is about twice as large as that determined by ga.o.- In other angular regions 
g gives less soft radiation than 5a.o.- The angular ordering approximation has the good 
feature that it gets the total amount of soft radiation right. 



This result follows from the original construction of Refs. [16]. We note, however, that 
the original construction involved only an integration over the azimuthal angle (j), while 
Eq. (7.12) requires an integral over both 6 and cf). We have also checked Eq. (7.12) by 
numerical integration. 

One should note that the theta function in g'a.o. restricts the emission angle of a soft 
gluon to be smaller than the angle between and ui, where A; is a parton that is color 



9a.o. 



(7.11) 




(7.12) 
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connected to parton /. If parton I is a quark, then there is only one choice for k. However, 
if parton I is a gluon, then there are two color connected partons. Then there are two 
contributions with separate angle restrictions. 

8. Conclusions 

In Ref. [1], we presented evolution equations that represent a leading order parton shower 
including quantum interference, spin, and color. We did not, however, present a way to 
implement the integrations implied by these equations in a fashion that would be practical 
for more than a few partons. The idea behind the evolution equations was to make just 
one approximation: that the virtualities in successive splittings are strongly ordered. 

Typical Monte Carlo event generators, such as Pythia [10], Ariadne [2], Herwig [11], 
and Sherpa [17], make additional approximations. In particular, they typically average 
over parton spins and take the leading term in an expansion in where A'^c = 3 is the 

number of colors. Our aim in this paper has been to work out how the general formalism 
could work as a practical calculation if we make the further approximations of averaging 
over parton spins^^ and of keeping only the leading order in We do, however, keep 

some aspects of quantum interference in that the interference graphs between the emission 
of a soft gluon from parton / and the emission of the soft gluon from another parton k are 
accounted for. 

The result is an algorithm that is similar to what is done in widely used parton shower 
event generators in that the calculation can be implemented as a Markov chain, as described 
in Sec. 6. The form of the evolution is perhaps most similar to that in the dipole showers 
of Refs. [4] and [5] and is also similar to the k± version of Pythia [3]. One can think of 
the basic object that splits as not one parton, but two partons, I and k, that are next to 
each other along a color string. This basic object is often referred to as a color dipole. 
When we incorporate the joint splitting of partons / and k, there is a contribution to the 
splitting probability that corresponds to the square of the amplitude for parton I to split. 
There is another contribution to the splitting probability that corresponds to the square 
of the amplitude for parton k to split. Then there are two contributions that correspond 
to the interference of these amplitudes. We reorganize the four terms into two terms. One 
is kinematically of the form for a splitting of parton while the other is kinematically of 
the form for a splitting of parton k. This is rather similar to the structure of the dipole 
subtraction scheme for next-to-leading order calculations proposed by Catani and Seymour 
[6], which has been implemented for parton showers in two recent papers [4, 5]. 

There are differences between the shower formulation used here and that in, say, the 
dipole showers of Refs. [4] and [5]. The splitting functions are different. In particular, 
we have separate formulations for the interference graphs (based on the simple eikonal 
approximation) and for the direct graphs, for which our splitting functions are quite directly 
read off from the Feynman graphs with a minimal approximation applied where an off-shell 
mother parton attaches to a hard scattering amplitude. The momentum mapping functions, 

^®More precisely, we average over the spins of a parton before it splits and sum over the spins of daughter 
partons. 
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which were presented in Ref. [1], are also different. They are similar to the Catani-Seymour 
momentum mappings in that they are systematically defined, invertible mappings, but they 
have the advantage that the form of the mapping depends on the parton index / but not 
on the index k of the partner parton. 

We have seen that the leading color, spin averaged shower of this paper has a structure 
similar to that implemented in standard parton shower event generators. In particular, this 
simple shower can be implemented using a Markov chain. The full shower formalism of 
Ref. [1] is more general than the simple shower in that parton spin and color correlations 
are included. We anticipate that the full formalism will be more difficult than the simple 
version to implement in a practical fashion. However, we anticipate that one can use the 
simple shower as a basis for a systematically improvable approximation to the full shower. 
The idea would be to start with the simple shower and provide parameters that remove the 
approximations gradually, so that the result is still approximate but the approximation is 
systematically improvable as computer resources allow. We expect to return to this subject 
in future papers. 
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A. The remaining splitting functions 

In this section we record the spin averaged splitting functions Wu for the cases in which 
fm+i 7^ g) which were not covered in the main body of the paper. We use the general 
definition (2.8) of Wu together with the formulas from Ref. [1] for the splitting amplitudes 

VI . 

We first consider a final state splitting with {/;, /m+i} = {g, Q, q} where g is a quark 
flavor and q is the corresponding antiflavor. A straightforward calculation gives 

u7 (sf -x ^ Stos ( 2 prD{pi,Q)-pm+i \ , . 

Wii{{f,p}.m+i) = TT— — ^2 1 + 7^——- ^2 • (^-^ 

[pi+Pm+lVy [pi+Pm+iY J 

For an initial state splitting with {fi, fi, fm.+i} = {g, Q, q}, we find 
TJr f[f -1 N Svras I , , f Pmi \^ 2 prn+i-D{puQ)-p.,n+i 



[Pl-Pm+lYy ^ ' \{pi - Pm.+l)-niJ [Pl-Pm+lY j 

(A.2) 

Here ni = p-Q for I = a and ni = px for / = b. The same result holds for an initial state 
splitting with {/;, fi, fm+i} = {g, q, q}- 

We consider next an initial state splitting with {fi, fi, fm+i} = {Q,§i,q}- A straight- 
forward calculation gives 

-nf-x \ fpi-ni {pi - pm+i)-ni Pm+i-D{pi,Q)-pm+i\ 

wii{{f,p}m+i) = — • (A.3) 

Pi-Pm+i \Pvni pi-ni Pl-Pm+l 
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Again, ria = pb and nb = PA- The same result holds for an initial state splitting with 

{fi,fi,fm+i} = {q,g,q}- _ 

This completes the analysis of Wu for cases in which fm+i / g- 
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